liS' 

!  r*  •■  :  y*  •- 

r'*  ' 

r  s.,j 


<0 

SO 


L  Iff-o  t?T 

FACULTY  DES  SCIENCES  DE  PARIS 


p:V  D‘.  - 


LA60RAT0IRE  DE  PHYSIQUE  DES  SOIIDES 


de  Recherche  Aasoci^e  au  C.N.R.S.  :  OPTIQUE  DES  SEMICONDUCTEURS 


NO  EOOAR -68-0016 

,f/?A  &/**  - 


SCIENTIFIC  REPORT  N°  1 


OPTICAL  STUDIES  OT  LATTICE  VIBRATION  IN  II- VI  SEMICONDUCTING  COMPOUNDS 


ifc^LbU  U  Lb 


THROUGH  THE  EUROPEAN  OFFICE  OF  AEROSPACE  RESEARCH,  OAR 


UNITED  STATES  ADR  FORCE  UNDER  GRAND  NUMBER  EO0AR-68-0016 


GRANT  N°  EOOAR  -68-0016 


SCIENTIFIC  REPORT  N°  1 


OPTICAL  STUDIES  OF  LATTICE  VIBRATION  IN  II -VI  SEMICONDUCTING  COMPOUNDS 


M.  BALKANSKI 


Laboratoire  de  Physique  des  Solides  de  la  Faculte  des  Sciences  de  Paris 
Equipe  de  Recherche  Associ6e  au  C.  N.  R.  S.  :  Optique  des  Semiconducteurs 


9,  Quai  Saint-Bernard  -  Paris  5feme 


This  research  has  been  sponsored  in  part  by  the 
AEROSPACE  RESEARCH  LABORATORIES 


Through  the  European  Office  of  Aerospace  Research,  OAR, 
United  States  Air  Force  under  grand  number  EOOAR -68-0016 
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IN  II-VI  SEMICONDUCTING  COMPOUNDS 


M.  Balk an ski 
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de  la  Faculty  des  Sciences  de  Paris 


1.  Introduction. 


In  recent  years  considerable  efforts  have  been  developped  to¬ 
ward  the  ’understanding  of  the  main  features  of  lattice  dynamics  of 
a  great  number  of  crystals  among  which  the  II-VI  semiconductor  com¬ 
pounds  have  a  good  position.  Neutron  scattering  technics  be^ing  in¬ 
applicable  to  most  of  these  compounds  the  experimental  determinations 
of  the  normal  mode  frequencies  are  essentially  done  by  infrared  spec¬ 
trometry  and  Raman  scattering.  Being  ionic  crystals  all  these  com¬ 
pounds  possess  first  order  dipole  moments  due  to  the  relative  dis¬ 
placement  of  the  constituting  ion3.  The  optical  vibrational  modes 
are  often  allowed  for  optical  transitions  in  view  of  specific  selec¬ 
tion  rules  and  give  the  possibility  of  direct  determination  of  the 
frequencies  of  these  modes  at  the  center  of  the  Brillouin  zone.  Those 
of  the  optical  modes  which  are  not  allowed  for  direct  optical  transi¬ 
tions  are  often  allowed  in  Raman  scattering  and  hence  available  far 
direct  determination  by  spectroscopic  technics. 

Normal  modes  frequencies  at  the  edge  of  the  Brillouin  zone  are 
determined  by  second  and  higher  order  processes  due  to  their  partici¬ 
pation  in  two  or  three  phonon  combinations  as  well  in  infrared  spec¬ 
troscopy  as  in  Raman  scattering. 

The  realistic  boundary  conditions  taking  into  consideration  the 
discontinuity  at  the  surface  of  the  crystal  bring  up  the  surface  modes 
of  vibrations  for  the  normal  lattice  and  the  localized  surface  modes 
for  a  crystal  with  foreign  atoms  absorbed  at  the  surface. 

Presence  of  imperfections  or  foreign  atoms  leads  to  the  appear¬ 
ance  of  localized  modes  of  vibrations  which  are  also  coupled  to  the 
radiation  field  by  means  of  resonant  absorption  in  the  infrared  region 
and  by  Raman  scattering. 
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Parallel  to  the  experimental  Investigations  theoretical  studies 
have  been  developed  and  the  dispersion  relations  calculated  for  same 
of  the  II-VT.  compounds.  The  theory  of  systems  with  high  impurity  con¬ 
centration  and  the  lattice  dynamics  of  the  alloys  is  now  under  way. 


2.  Crystal  Structures. 

The  I I- VI  semiconducting  compounds  crystallize  in  a  variety  of 
polymorphic  modifications,  the  principal  structure  types  being  : 

a)  Cubic  zinc-blende  (sphalerite)  structure  \sed  on  the  cubic  space 
group  T^  s 

ZnS,  ZnSe,  ZnTe,  CdTe,  HgS,  HgSe,  HgTe  crystallize  preferen¬ 
tially  in  the  zinc-blende  structure. 

Each  unit  cell  contains  2  atoms  and  the  crystal  will  therefore 
possess  6  phonon  branches. 

b)  Hexagonal  wurtzite  (zincite)  structure  based  on  the  hexagonal 

4 

space  group  Cgy  s 

ZnO,  ZnS,  CdS,  CdSe  crystallize  preferentially  in  the 
wurtzite  structure. 

Each  unit  cell  contains  4  atoms  and  the  crystal  will  therefore 
possess  12  phonon  branches. 

A  schematic  representation  of  these  two  crystals  crystals  struc¬ 
tures  is  given  in  Pig. 1 . 


Figure  1.  Schematic  representation  of  Blende  and  Wurtslte  crystal  structures. 


On  the  left  side  is  shown  the  cubic  structure  where  one  can 


distinguish  three  different  successions  of  layers  A,  B,  and  C,  with 
alternating  atoms.  On  the  right  side  is  the  hexagonal  structure  where 
one  finds  only  two  different  successions  of  layers  A  and  B  respectively. 

Or.  Pig.  2  are  schematized  the  two  elementary  cell3  containing 
two  atoms  for  the  cubic  structure  shown  on  the  left  and  four  atoms 
for  the  hexagonal  structure  on  the  right. 


Figure  2.  Representation  of  elementary  cells  : 

a)  left  cubic  structure 

b )  right  hexagonal  structure 


The  first  nearest  neighbours  in  both  structures  are  at  the  cor¬ 
ners  of  a  regular  tetrahedron  related  to  those  of  the  group  IV  semi¬ 
conductors.  For  this  reason  many  physical  properties  are  analogous 
in  both  types  of  crystal  lattices. 
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The  two  respective  environments  in  the  cubic  and  hexagonal 
structure  are  shown  in  Pig. 3. 


Figure  3.  Environments  In  the  cubic  left  end  hexegonel  right  structure. 

The  hounding  in  the  crystals  of  the  II-VT  compounds  spares  the 
complete  gap  from  ionic  to  covalent. 

The  first  Brillouin  zone  deduced  from  this  structure  is  shown 
in  Fig. 4  respectively  for  the  blende  and  the  wurtzite. 


Figure  4.  First  Brillouin  sons  for  Blende  and  Wurtslte. 
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3.  Dispsrsion  Relationo. 


In  ths  harmonic  approximation  where  the  potential  energy  of  a 
cryatal  ia  expanded  in  power  of  the  diaplacament  of  the  atoms  from 
their  equilibrium  positions  and  the  axpanaion  ia  broken  off  the  quad¬ 
ratic  terms*  the  equation  of  motion  of  the  atom  k  with  maos  in  the 
unit  cell  1  subject  to  ths  diaplacament  uq  with  a  =  x,  y,  z  is  given 


by 


*  (n'i  u  I1')  = 
*a^kk’'  -Vk,; 


(3.1) 


where  #  ie  the  force  constant  matrix  dspending  only  on  the  distance 
between  cells  j(l)-»(l')  due  to  the  periodic  symmetry  of  the  lattice. 


The  normal  modes  and  thair  frequencies  arc*  solutions  of  ths  set 
of  equations  for  all  unit  calls  of  the  crystal.  Because  of  ths  period¬ 
ic  symmetry  the  modes  have  wave  like  phase  variation  from  cell  to  call 
and  specific  relative  motion  in  the  cell. 


Ths  dynamical  matrix  being  of  the  farm 


( 3.2) 


ths  system  of  3n  linear  and  homogenious  equations  has  non-vanishing 
solutions  only  if  the  determinant 


CaP^k  k’  ^  "  w  6 


kk'  ap 


=  0 


(3.3) 


vanishes. 


This  "characteristic"  equation  gives  ths  chractsristic  frequen¬ 
cies  wCj)  whose  graphic  representations  are  ths  dispsrsion  curvss  of 
the  crystal. 


A  certain  number  of  investigations  have  bsen  developed  aiming 
the  knowledge  of  the  dispersion  curvss  of  ths  II-VI  compounds.  Some 
of  them  ’  bassd  on  the  shsll  modsl  develop  ths  formal  discussion 
without  actually  doing  the  calculations.  Calculating  the  dispsrsion 
relations  of  blends  or  wurtzits  structures  in  the  shell  model  would 
involve  a  large  number  of  parameters  which  make  such  calculations 
uaaninglsss. 

3 

Marten  has  developed  a  series  of  studies  in  which  the  lattice 
vibrations  of  blends  and  wurtzits  structures  are  investigated  an 
microscopic  basis  considering  ths  dynamical  matrix  without  and  with 
long  rangs  Coulomb  interactions.  Here  again  the  numerical  calcula¬ 
tions  havs  not  been  carried  out  throughout  the  whole  Brillouin  zone 
and  only  the  optical  modes  at  the  center  of  the  Brillouin  zone  have 
been  deduced. 
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Recently  Nusimovici  and  Birman  ^  have  performed  the  explicite 
calculations  of  the  dispersion  relations  for  CdS-wurtzite  structure. 
The  single  phonon  frequencies  have  been  used  to  calculate  the  farce 
constants.  The  model  considers  short  range  interactions  correspond¬ 
ing  to  central  forces  between  first,  second  and  third  neighbours  and 
to  central  angular  forces  between  first  and  second. neighbours,  as 
well  as  long  range  Coulomb  interactions  depending  on  a  parameter 
having  the  dimensions  of  a  spring  constant. 

The  wurtzite  is  supposed  to  be  ideal  and  from  the  point,  of 
view  of  the  short  range  interactions  the  anisotropy  is  represented 
by  the  third  neighbour  forces  as  shown  on  Pig. 5. 


Figure  5.  Anisotropy  deduced  from  short  range  Interaction  In  Wurttlta. 
The  short  range  potential  is  written  as 


where  X,  p,  v,  and  6,  correspond  to  central  forces  between  first, 
second  and  third  neighbours,  kg,  kg,,  krg  end  krg,  correspond  to 
non-central  angular  forces  between  first  and  second  neighbours. 
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PHONON  ENERGY 
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A  calculation  done  for  an  hypothetical  cubic  CdS  shows  the  dif¬ 
ferences  between  the  two  phases  :  wurtzite  and  blende. 


Figure  7.  Dispersion  relations  for  Blende  for  the  M,  r  and  L  points 

of  the  Brlllouln  Zone. 

An  important  check  of  their  calculation  is  given  by  the  two- 

phonon  optical  data  obtained  by  infrared  and  Raman  spectra.  Table  I 

compares  the  two- phonon  Raman  scattering  in  CdS  calculated  by  Nusi- 

5 

movici  and  Birman  with  the  ones  obtained  by  Tell,  Damen  and  Porto  . 

In  more  recent  work  ^  Nusimovici  and  Birman  made  a  direct 

7 

comparison  between  infrared  spectra  obtained  experimentally  and 
the  calculated  two- phonon  density  of  modes. 
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Table  I 


Two- phonon  Raman  scattering  in  CdS-wurtzite 

Experiment 
{ cnf 1 ) 

Interpretation 

Calculation 
( cm" 1 ) 

97 

VO 

P-. 

CM 

90 

2M4 

86 

207 

© 

"3 

201 

M2 

9 

M3 

202 

K2 

© 

K 

209 

*5 

9 

*3 

210 

328 

“ 

© 

« 

323 

*3 

© 

*4 

328 

© 

331 

347 

A1 

9 

342 

«, 

© 

*2 

345 

«, 

9 

344 

M2 

© 

M2 

346 

«, 

© 

», 

342 

X, 

© 

M3 

346 

K2 

© 

347 

K2 

© 

K2 

346 

*3 

© 

*3 

544 

K3 

© 

K2 

346 

364 

Ai 

© 

363 

K 

© 

K2 

368 

K 

© 

*3 

369 

Mi 

© 

“l 

368 

556 

2K 

558 

2M3 

550 

605 

2ri 

596 

2T5(D 

616 

4.  Photon- Phonon  Interaction. 


One  of  the  main  possibility  to  obtain  experimental  precisions 
on  the  diepereion  relations  and  particularly  the  normal  modee  fre¬ 
quencies  at  the  center  and  the  edges  of  the  Brillouin  zone  are  infra¬ 
red  spectroscopy  and  Raman  scattering.  For  the  II-VI  compounds  of 
blende  and  wurtzite  structures  these  two  methods  are  complementary 
because  of  particular  selection  rules  based  on  symmetry  considera¬ 
tions. 
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The  optical  absorption  is  given  by  the  transition  probability 
proportional  to  a  matrix  element  which  is  different  from  zero  only 
when  energy  and  momentum  conservation  rules  are  obeyed.  The  matrix 
element  for  optical  absorption  depends  on  the  dipole  moment  operator 
Y-  It  is  then  easy  to  deduce  by  group  theory  arguments  the  optically 
active  normal  modes  in  blende  and  wurtzite  structures.  These  are  of¬ 
ten  different  from  that  active  in  Raman  scattering  while  in  those 
cases  the  matrix  element  is  related  to  the  polarizability  tensor  H. 

4.1.  Infrared  Dispersion. 

When  the  electromagnetic  field  acts  on  the  crystal  there  is  a 
coupling  between  the  elastic  waves  giving  the  phonon  fields  and  the 
electromagnetic  waves  giving  the  photon  fields. 


The  electromagnetic  field  is  given  by  Maxwell’s  equations  : 


CSJ  8  =  “  £  , 

Mi  Emt =  -  i  i  • 

div  2=0,  (4.1) 


div  a  =  0  , 


and  fi  *  fijnt  +  4Ttp  .  (4.2) 

The  phonon  field  is  described  by  the  induced  polarization  duo 
to  the  displacement  in  opposite  directions  of  the  two  sublattices 
farmed  by  charged  ions  and  the  internal  electric  field  Jg  equal 
to  the  macroscopic  electric  field  in  the  crystal, 


£  =  H  e„  U  + 


o  -1 
00 


4K 


*int 


(4.3) 


and  also  by  the  equation  of  motion 

N  0  q 

M  a  +  (Mm02 - -  — ) U~3q,  Sint  •  (4.4) 

The  solutions  of  the  coupled  system  of  equations  representing 
the  radiation  field  and  the  phonon  field  ,-ive  the  infrared  dioporsion 
for  a  given  crystal. 


A  detailed  analysis  of  infrared  dispersion  for  blende  and  wurtz- 

8 

ite  structures  has  been  carried  out  by  Le  Toullec  ;  this  chapiter 
summarize  the  results. 


11 


4.1.1.  Blende. 


In  a  rigid  ion  approximation  when  the  ionic  charge  is  centered 
on  the  ion|  the  dynamical  matrix  can  be  written  in  two  terms  : 


C  i  4 
Co0lk  k 


\  -  CN  ,  3  )  +  cc  (  S  ) 
' '  a{rk  k' '  0aPlk  k' ’ 


(4.5) 


The  first  term  concerns  the  short  range  interaction  and  the  second 
the  Coulomb  interaction  due  to  long  range  forces.  The  polar  modes 
which  are  infrared  active  arise  from  these  Coulomb  forces. 


The  first  term  gives  all  the  nrrmal  modes.  For  q  =  0.  the  secu- 

p 

lar  equation  (C-orl)  gives  two  distinct  roots  threefold  degenerate  : 


to  =0 

ac 


rmi+m2  y—  N  /  1  \-i 


1/2 


(4.6) 


The  displacement  vectors  corresponding  to  to  are  y  =  y  ,  with  the 

&C  1  c 

two  sublattices  moving  i.i  the  same  direction,  those  corresponding  to 
&>op  are  y  1  =  )y^  >  with  the  two  sublatticea  moving  in  opposite 

directions.  To  the  optical  modes  corresponds  on  induced  dipole  moment 
equal  to 


M  =  e#  y 


Vm2 

m„ 


(4.7) 


In  a  rigid  non-polarizable  ions  approximation  the  Coulomb  potential 
energy  can  be  written 


•  -■Jr 


e«  e,  * 
k  x1 


2  ll7  r(J)  -  p(J|) 
kk'  *  * 


(4.8) 


and  the  coefficients  C  deduced  from  ®  for  q  =  0  are 


C 


C  (  0 
apvk  k' 


)  =  - 


4ft  N  , 


(4.9) 


whore  N  is  the  number  of  ion  paires  per  unit  volume.  If  we  hav',  the 
propagation  vector  y  along  the  Oz  axes  we  obtain  now  two  optical 
branches  instead  of  one.  The  first  is  twofold  degenerate  cor¬ 
responding  to  vibrations  perpendicular  to  y  : 


411  N  e^ 

3M 


where 


al 


(4.10) 


12 


where  u)q  is  the  optical  frequency  when  the  Coulomb  farces  are  neglect¬ 
ed  and  the  second  non-dsgenerate  corresponding  to  vibrations 
parallel  to  a  : 


871  N  e  2  471  N  e  2 

“lO  =  “o  +  3M  =  “to  *  5 


(4.11) 


Comparjd  to  equation  (4.4)  we  see  that  these  two  equations  differ  by 

2 

the  effective  charge  a  instead  of  e  e_  due  to  the  fact  that  we.  have 

1  8  1 

considered  rigid  non-polari zable- ions. 


The  displacement  in  the  case  of  polarizable  atoms  will  be 


T  «int 
M(«To-u>2) 


(4.12) 


Substituting  in  equation  (4.3)  and  subsequently  in  equation  (4.2)  we 
obtain  the  equation  of  optical  dispersion  of  the  crystal 


e(a)  =  e 


471  N 

M(o.t2-(02) 


(4.13) 


and  in  the  limit  u>  =  0  we  have  the  static  dielectric  constant 


471  N  e  ‘ 

C°s£"+^r 

where  from 

,  ,  (VeJ  ^0 

e(u.)  =  e#Q  +  - 5 - — 

“TO"  u 

two  cases  can  be  distinguished. 

First  case  :  e  =  0  ;  y  // a  // £in^.  '•  longitudinal  wave. 

(eo-E«^T0 
e  =  0  =  e  +  — -a - Sr 

00  d  d 

-  0) 

TO 


hence 


2  2  2  471  N  e  c 

“L  "  UT0  ^  "  “TO  +  M  e~ 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


This  frequency  is  independent  of  the  wave  vector  q  and  is  equal  to 
the  longitudinal  optical  frequency  of  the  system  without  the  external 
field  being  applied. 

Second  case  :  J4  1  a  1  £•  t  !  transverse  wavs. 


From  Maxwell's  equations  one  deduces 


.2  "  Eoo 


(Ve>T0  ,  * 

+  - 2 - 2~  = 


“to  -  “ 


(4.18) 


The  incident  light  interacts  strongly  with  transverse  optical 
phonons  of  the  crystal.  The  normal  mode  of  the  coupled  system  ; 
photon  +  phonons  (TO)  is  called  polariton.  The  dispersion  relation 
of  polaritons  describes  the  phenomena  of  optica1  dispersion.  The 
dispersion  u  =  f(q)  curve  could  be  easily  deduced  if  one  considers 
the  Lyddane-Sachs-Teller  relation 


“lo  "  “to  (ec/eJ 


Equation  (4.18)  is  now  written 


d  =  ~  eoo  “ 


“2~  “LO 

“2  -  “to 


(4.19) 


(4.20) 


The  dispersion  curve  of  polariton  gives  w  =  f(q)  for  a  system  of 
coupled  photons  plus  optical  phonons  is  represented  in  Pig. 8. 


Figure  8.  Dispersion  curve  of  optical  phonons  coupled  with  photons. 
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4. 1.1. a)  Introduction  of  Damping. 

In  the  harmonic  approximation  all  phonons  are  independent  from 
each  other  and  therefore  the  equations'  of  motion  are  linear.  In  real¬ 
ity  the  system  of  equations  is  coupled  through  the  unharmonic  terms 
of  the  potential.  Part  of  the  electromagnetic  energy  brought  into 
the  TO  by  the  incident  radiation  field  is  transferred  to  other  pho¬ 
nons.  The  introduction  of  a  damping  constant  in  the  phenomenological 
equations  expresses  this  fact  and  the  equation  of  motion  becomes 


“  1»  +  m  Y  U  +  m  u>T0  jj  =  erj,  gint 


Equation  (4.15)  should  now  be  written 

,  \  ^ Eo~E«' ^TO  ,  i2 

e(o))  =  £„  +  — p - 5 -  =  •  +  liC^  * 

WT0  ‘  “  "  iy<Ji 


(4.21) 


(4.22) 


The  real  and  imaginary  parts  of  the  dielectric  constant  are  then  : 

2  .  2\  .  2 


r  _  2,2  (yO^TQ-^Ho 
f~  "  00  (u>t?-“>2)2+  y2«2 

The  frequency  pole  of  e(u)),  can  be  written 


(4.23) 


(4.24) 


U)m  =  0)  +  10), 

T  r  1 


with 


0) 


1/2 


TO 


and 


TO 


(4.25) 


(4.26) 


w  =  -  1  . 

i  2 


The  real  part  of  0)^,  represents  the  resonance  frequency  and  dif¬ 
fers  very  little  from  uT0  since  (y/u>tq)  1 

The  imaginary  part  of  measures  the  damping. 


“t 


is  a  root  of  e(u))  =  0  and  can  be  written 


WL  =  Culo’^2/4)]  -  [iT/2] 


1/2 


(4.27) 


L 


The  infrared  dispersion  curves  for  the  transverse  modes  are 
rep:  esented  on  Pig. 9. 


u 


Figure  9.  tnhariaonic  diaper aion  curve  for  the  TO  mode. 


Figure  10.  Anharoonic  optic  conatanta. 


_ C,=  l>2-k2 

_  ti,z  2n  k 
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Figure  11 .  Anharaonic  dielectric  Constanta 
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4. 1. 1.b)  Experimental  results  from  infrared  spectro¬ 
scopy. 

Dispersion  in  blende  type  crystals  has  been  studied  by  many 
authors.  Ih  recent  years  systematical  studies  have  been  developed 
as  well  by  reflectivity  as  by  absorption  measurements  in  the  infra¬ 
red  region.  Receizt  results  for  blende  are  summarized  on  Table  II. 

Table  II  -  Blende 

Frequencies  of  the  optically  active  modes  in 
Blende  type  crystals  determined  by  infrared 
spectroscopy. 


WT0 

UL0 

eo 

£ 

00 

References 

ZnS 

274 

350 

8.3 

5,07 

a,  b 

ZnSe 

207 

242 

7,6 

5.4 

c 

209.5 

250 
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5.75 

d 

209 

257 

9.2 

6.1 

e 

9.1 

f 

8.66 

5 

5.90 

h 

ZnTe 

177 

205 

9.1 

6.7 

c 

10. 1 

f 

9.67 

g 

7.28 
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CdTe 
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4.1.2.  WUrtzite. 

4. 1.2. a)  Short  ranae  forces. 

The  unit  cell  contains  4  atoms  2  anions  and  2  cations,  12  branches 
are  therefore  expected  in  the  dispersion  curves. 

The  following  relation  between  C  coefficients  results  from  the 
crystal  symmetry  i 


CS  (  0  )  -  r$  (  G  W  CN  (  o  ) 
xx^k  k* '  Wk  k*'  +  Gzzlk  kJ  * 


(4.28) 


One  obtains  6  different  solutions  of  the  secular  equation.  One  is 
threefold  degenerated  and  corresponds  to  the  acoustical  branches. 

Prom  the  5  other  solutions  for  the  optical  frequencies  2  are  ncn-degensn- 
ate,2  are  twofold  degenerate  and  one  threefold  degenerate  due  to  an 
accidental  degeneracy  for  • 

The  relative  motion  of  the  four  atoms  constituting  the  unit 
cell  is  schematically  re;  esented  on  Pig. 12. 


41  j 


t°inO*.sro  (E,.r.)  ujli 


J  W2II  (B, ,  !■;>  (E„0  C021 


W3,  .  ri)  fcO  "31 


figure  12.  Relative  motion  of  the  four  atone  In  the  hexagonal  unit  cell. 
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In  the  optical  modes  of  vibration  for  the  frequency  the  sub¬ 
lattices  with  the  same  kind  of  atoms  (l)  and  (3)  are  moving  in  the 
same  direction  and  the  sublattices  with  atoms  (2)  and  (4)  in  opposite 
directions’.  The  displacement  vectors  are  related  by 


u  (l)  =' - -  u  (2)  =  u  (3)  - - —  u  (4) 

~a  m1  -a  -a  m^  -,av  ' 


(4.29) 


This  kind  of  displacement  results  into  a  total  dipole  moment 
per  unit  cell  s 

m  +m 

M=2e*»1  *  (4*30) 


For  the  nodes  of  frequencies  Ug  and  the  displacement 
vectors  are  of  the  same  type  but  the  induced  dipole  moment  per  unit 
cell  is  always  equal  to  zero  and  therefore  these  modes  are  not  optic¬ 
ally  active. 

Fig. 13  shows  the  relative  position  of  the  normal  mode  frequen¬ 
cies  near  the  center  of  the  Brillouin  zone  deduced  only  from  the 
consideration  of  the  nearest  neighbours  short  range  interaction. 


“2J. 

1 


*VlWU 


Figure  13.  Position  of  the  normal  mod*  frequencies  near  q 
of  short  range  Interaction  only. 


■  0  in  the  case 
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4 . 1 . 2 . b )  Long  range  Coulomb  forces . 

Q 

The  coefficient  C  deduced  from  the  Coulomb  potential  energy- 
can  be  written  in  the  case  of  wurtzite  : 


„C  ,  0  v  N  ek  V  / _kk 1  •  411  \  V  ,.s 

“Sik'  (\v>^  '  “r  i2  ’ 

kk* 

are  Lorentz  factors  in  anisotropic  crystals  and  their  numerical 
values  are  known. 


Introducing  Coulomb  interaction  into  the  determination  of  the 
normal  mode  frequencies  one  have  also  to  distinguish  between  the 
case  where  the  propagation  vector  q  is  parallel  or  perpendicular  to 
the  optical  axis  q.  The  modes  corresponding  to  atom  vibrations  per¬ 
pendicular  to  the  propagation  vector  labelled  TO  and  these  to  parallel 
vibration  LO,  //  and  1  refers  to  the  vibration  of  the  atoms  with  re¬ 
gard  to  the  q  axis.  The  optically  active  triply  degenerate  mode 
without  Coulomb  interaction,  splits  now  into  the  following  frequencies 


For  aifi  * 

“LO  =  “1  +  8,47  T 


“toj_=  “I 


-  4.09 


Ne£ 

if 


1(0v/ 


For  q  //q  i 

Wg  ^  Wg  ^ 

•wy-  “V/+  8-17ir  5  “soi uil  "  4,09  T  * 

Here  N  is  the  number  of  unit  cell  per  cm^,  e  is  the  electronic 
charge,  and  M  the  reduced  mass  :  M  =  n^m^/Cm^m^)  .  These  results  are 
schematically  represented  in  Fig. 14. 


4.l.2.c)  Infrared  dispersion  curves. 

Wurtzite  being  on  uniaxial  crystal,  because  of  the  anisotropy, 
the  infrared  dispersion  curves  depend  on  the  direction  of  propagation. 
The  solutions  of  the  coupled  equations  giving  the  radiation  field  and 
the  phonon  field  are  now  of  the  form  t 

(n2-C|)  [(n2-e^)(n2-€^)-n2sin2e(n2-e^)-n2cos2e(n2-e|)]  »  0  .  (4.32) 

0  is  the  angle  between  the  propagation  vector  and  the  optical  axis. 
Two  types  of  solutions  represented  in  Fig. 15  appears. 
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Figure  14. 


with 

coulo»b 

forces 


without 

couloab 

forces 


with 

coulowb 

forces 


U10II 


“toi 

<*1211 


"31 


CJJA 


Splitting  of  dispersion  curves  near  q  *  0  when  Coulomb 

forces  ere  taking  Into  account. 


a)  on  the  left  the  propagation  vector  Is  perpendicular  to  the  c  axle 

b)  on  the  right  the  propagation  vector  le  parallel  to  tha  c  axle. 


ordinary  extraordinary 


Figure  15.  Dispersion  curve  dependence  on  the  direction  of  propagation 
•)  //  e 

b)  q  1  c 
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A)  Ordinary  solution. 


n  -  Ei  =  0  . 


This  solution  is  indeperdent  of  the  angle  6  (Pig. 4a).  The  cor¬ 
responding  dispersion  rslation  ls 


1 


2  "Hal 


(4.34) 


B)  Extraordinary  solution. 


This  second  solution  corresponding  to  the  second  factor  of  the 
equation  (4.32)  depending  on  the  angle  0.  Two  extreme  cases  can  bs 
considered  : 


0  =  0  i.e.  a//  £  and  0  =  ■?  i.e.  a  1 


-  First  case  :  a  //  £  5  9  *  0  . 

The  equation  (4.32)  becomes  now  : 

(4.35) 

z  0  gives  u=  »  solution,  independent  of  a* 

The  infrared  active  longitudinal  phonon  do  not  influence  the 


=  0 


infrared  dispersion  in  the  considered  direction. 

.2 


n  -E|  =  0  gives  the  second  transverse  solution  (Fig. 4b)  with  the 
same  dispersion  rslation  as  equation  (4.34). 

-  Second  ca3S  :  a  J.  £  5  9  =  ■?  . 


Equation  (4.32)  becomes 


ej^(n-e^)  =  0 


(4.36) 


where  =  0  with  u  * 


n  -Ey  =  0  with.ths  dispersion  solution 

2  _ _ 

..2  .2 

V/ 


2  1  2 
£  ?  U 
cz  “  o-ufOj 


(4.37)  (Fig. 4b) 


ffe  have  see’i  that  urjoj^  is  a  transverse  optical  mode  £  1  a  with 
£  1  £  >  ,dTQy  a  transverse  optical  mods  £  J.  a  with  £  /  £  .  The 
solutions  of  equation  (4.32)  axe  given  on  the  Pig. 4c. 


Experimental  determination  of  these  frequencies  are  made  by 
the  following  two  methods. 


-  ***** 
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The  reflectivity  is  measured  under  weak  incidence  angle  (0~i5°). 
The  curves  representing  the  refractive  index  n(w)  and  the  extension 
coefficient  K (u)  are  calculated  by  means  of  Kramer s-Krbnig  inversion 
formula. 


The  frequency  is  the  frequency  for  which  the  quantity 
(2nKu)  has  its  maximum. 

The  frequency  is  the  frequency  for  which  one  has  the  rela¬ 
tion  n2-K^  =  0  with  (&i/bu>)  ^  0  . 


In  order  to  be  able  to  determine  the  four  frequencies  a^Q 
“TOJ^  »  **1*0//  •  sufficient  to  prepare  a  sample  with  g 

axis  perpendicular  to  the  plane  of  incidence  and  to  use  a  polarized 
light  beam  either  parallel  to  the  g  axis  in  which  case  one  obtains 
oiLO^and  wto^i  or  perpendicular  to  the  g  axis  in  which  case  one  ob¬ 
tains  811(1  «i0J.  * 


2)  Thin  films. 

When  transmission  or  reflectivity  measurements  are  performed 
on  thin  films  under  oblic  incidence  and  with  a  light  beam  polarized 
in  the  plane  of  incidence,  the  transverse  optical  as  well  as  the 
longitudinal  optical  modes  interact  with  the  infrared  radiation  and 
hence  their  frequencies  can  be  measured  directly.  In  the  case  of  an 
uniwTfal  crystal  if  the  thin  film  is  such  as  to  have  the  g  axis  per¬ 
pendicular  tothe  substrate  one  measures  then  0)*  and  the  <o<poj.  fre¬ 
quencies.  This  is  the  case  of  CdS  thin  films  evaporated  on  silicon 
single  crystal  substrates. 

If  the  g  axis  is  parallel  to  the  surface  of  the  layer  it  is 
then  possible  to  observe  wlqi • 


The  different  experimental  situations  in  which  measurements 
have  been  made  by  Le  Toullec  ®  are  listed  in  Pig. 16. 


ura  16. 


Comparison  of  the  optical  frequenciee  of  bulk  aamplee  and 
thin  filma,  under  oblique  incidence  and  with  polarised  light 
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The  experimental  results  are  summarized  in  Table  III. 
Table  III  -  Wurtzite. 

Frequency  of  the  optically  active  modes  in 
wurtzite  structure  crystals  determined  by 
infrared  spectroscopy. 


“TO ^ 

cm- 1 

“TOI 
cm- ' 

“LO/ 

cur1 

“LOI 

cm- ' 

E°// 

E°1 

*"// 

EooJ^ 

References 

232 

240 

298 

302 

8.9 

8.40 

5.38 

5.3 

a 

233.6 

241.2 

299 

305 

8.9 

8.50 

5.38 

5.29 

b 

CdS 

9.53 

9.02 

c 

235 

242 

296 

303 

6.42 

8.37 

5.32 

5.32 

d 

167 

171 

209 

211.5 

9.95 

9.1 

6.05 

5.95 

a 

166 

172 

211 

210 

10. 16 

9.29 

6.3 

6.2 

e 

CdSe 

6.05 

5.95 

f 

10.20 

9.53 

c 

168 

171 

211 

211 

9.6 

9.3 

6.1 

6.1 

d 

References 

a)  Le  Toullec,  Thesis,  Paris. 

b)  A.  Mitsuishi,  Paper  presented  at  "U.S. -Japan  Cooperative  Seminar 
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(1963). 

d)  H.W.  Verleur  and  A.S.  Barker  Jr.,  Phys.  Rev,  155 .  750  (1967). 
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f)  Shiozawa  and  al.,  Research  on  II-VI  compound  semi conduct ore,  6  11 
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Aeronautical  Research  Laboratories,  Wright-Patterecn  Air  Force 
Base,  Ohio. 


4.2.  Light  scattering. 

The  inelastic  ecattering  of  the  phonons  by  the  crystal  is 
generally  referred  to  ae  Raman  ecattering  if  optical  phonone  are 
involved  in  the  processus  and  as  Brillouin  scattering  if  the  acous¬ 
tical  branches  are  involved.  Phonone  may  be  absorbed  or  emitted  and 
accordingly  stokes  or  antistokee  linee  will  be  observed  on  the  long 
wave  eide  or  the  short  wave  side  of  the  incident  radiation.  The  polar- 
ieibility  of  the  cryetal  ie  composed  of  two  parte  ;  one  due  to  the 
nuclear  displacements  and  the  other  to  the  consequent  electronic 
polarization. 


25 


¥ 

£~“ 


In  first  order  the  diffused  intensity  for  the  stokes  Raman  line 
in  the  harmonic  approximation  is 


stokes 


=  K  v 


diff 


nTV+1  <  . 

t«a 


4lt!) 


n  is  the  number  of  phonons  occupying  the  initial  p  state,  is  the 

P  i  th  ° 

oscillator’s  frequency,  e  is  the  a  component  of  the  unit  polariza- 

O  d iff  th 

tian  vector  of  the  incident  radiation,  e“  is  the  b  n  component  of 

the  polarization  unit  vector  of  the  diffused  light,  -  ^“ab^^e 

is  the  variation  of  the  component  of  the  polarization  tensor  for 

the  vibration  with  normal  coordinates  q  . 

e 

Symmetry  considerations  and  group  theory  arguments  define  what 
are  for  each  normal  mode  the  polarizibility  tensor  components  which 
have  changed  during  the  vibration.  It  can  be  easily  shown  that  in  the 
blende  type  crystals  the  two  mcdes  at  the  canter  of  the  Brillouin  zone 
are  Raman  actives. 


Experimental  results  giving  the  frequencies  of  the  Raman  active 
optica  as  in  blende  structure  are  listed  in  Table  IV. 

Table  IV 

Observed  frequencies  of  Raman  active  optical  modes 
in  blende  structure  crystals. 


“to 

“lo 

References 

ZnS 

274 

349 

a)  L.  Couture-Mathieu,  J.P.  Mathieu,  C.R.  Acad.  Sci. 
Paris  22£,  371  ( 1953). 

205 

253 

b)  S.S.  Ultra,  J.  Phys.  Soc.  Japan  Suppl.£1,  61  (1966) 

ZnSe 

o 

CM 

251 

c)  U.  Krauzman,  C.R.  Acad.  Sci.  Paris  264, 1117  (1967) 

205 

250 

d)  W.  Taylor,  Phys.  Letters  556  (1967). 

ZnTe 

177 

206 

c)  M.  Krauzman,  C.R. Acad. Sci.  Paris  264.  1117  (1967) 

174 

203 

d)  W.  Taylor,  Phys.  Letters  556  (1967>. 

4.2.2.  Wurtzite. 


Selection  rules  for  Raman  scattering  in  wurtzite  structure  crys- 

q 

tals  have  also  been  deduced  by  group  theory  arguments  .  The  Raman 
active  modes  in  wurtzite  are  the  dollowing  : 


lOfjpS, )  |  TO^U,)  }  T02(E2)  ;  T03(E2)  5  LO^E,)  ;  LO^U,) 
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In  Table  V  are  given  the  experimental  results  determining  the 
frequencies  of  the  Raman  active  modes  in  ZnO  and  CdS  single  crystals. 


Table  V 

Observed  frequencies  for  Raman  active  optical  modes 
in  wurtzite  structure  crystals. 


“LO^ 

“LOy 

wT02 

“TO^ 

“«>1/ 

^03 

References 

ZnO  583 

574 

437 

407 

380 

1 01 

a) 

538 

588 

333 

438 

420 

180 

b) 

CdS 

325 

212 

256 

85 

c) 

305 

305 

252 

235 

228 

44 

d) 
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b)  S.S.  Mitra  and  J.I.  Bryant,  Bull.  Am.  Phys.  Soc.  JO,  333  (1965). 

c)  H.  Poulet  and  J.P.  Mathieu,  Ann.  Phys.  (Paris)  549  (1964). 

d)  B.  Tell,  T.C.  Damen  and  S.P.S.  Porto,  Phys.  Rev.  144.  771  (  1966 ) . 


Multipnonon  Interaction. 

The  normal  modes  frequencies  at  the  edges  of  the  Brilluoin  zone 
can  be. deduced  from  infrared  absorption  due  to  multiphonon  interaction. 
Generally  all  phonons  can  participate  to  multiple  processes  ;  the  ac¬ 
tual  transitions  obey  selection  rules  deduced  from  symmetry  considera¬ 
tions. 

The  theory  of  the  interaction  of  the  radiation  field  with  two 

or  more  normal  vibration  modes  has  been  given  in  terms  of  second  or 

higher  order  dipole  moments  or  as  due  to  the  presence  of  aiharmonic 

terms  in  the  potential  energy  associated  with  lattice  vibrations  !l 

'2 

and  also  as  a  combination  of  both 


In  the  case  of  optical  absorption  due  to  double-phonon  processes, 

the  contribution  to  the  imaginary  part  of  the  dielectric  constant  due 

to  anharmonic  terms  higher  than  second  order  in  the  potential  energy 
13 

can  be  written  : 


e"(«)  =  rcK  r 
33' 


2  IT 

n,r 

2  2 

w) 

TO"U33’ 

Vi' 

\  *  (n^  +  j) 


|  grad  d>(a)| 


ld=W 


33' 


+  higher  order  terms  , 


Here  =  e*(N/M)-l/2  ,  bj;j,T  = 


^Qj(q)  &Q3,(q')  ^Qto 
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is  the  coupling  fo re  between  the  TO  mode  and  the  phonons  of  wave 
vector  <j.  *-■  ±  j. 


,  n^  is  the  thermal  mean  value 


M  ~  ~  *y  '  V  =  “3  *  “3*  , 

of  the  quantum  number  n^  .  The  factors  (n^.+n^+l)  and  (n^,-n^)  give 
the  temperature  dependence  of  the  additive  and  subtractive  absorption 
bands. 


The  summation  is  taken  over  the  surface  w(a)  in  the  reduced 
Brillouin  zone  for  which  u(a)  =  critical  points  for 

this  "combined"  density  of  states  give  rise  to  bands  in  the  absorp¬ 
tion  spectrum. 

The  contribution  to  the  imaginary  part  of  the  dielectric  con¬ 
stant  due  to  second  or  higher  dipole  moments  can  be  written  : 


e"(«)  =  n2#  JZ 


*(n> 


ni»  +  2  ~  ^n.i *  2L 
I  grad  “(a)  I  ^ 


da 


33' 


+  higher  order  terms. 


Here  Jf(^  ”*, )  =  b2H/bQ(^) .bQ(“^)  is  the  change  of  effective  charge 
associated  with  a  given  phonon  (^)  due  to  a  second  phonon 

Thus  independently  of  the  nature  of  the  processes  with  which 
one  is  concerned  the  characteristic  variation  of  the  absorption  coef¬ 
ficient  as  a  function  of  the  temporature  is  the  same. 


.Analysis  of  the  absorption  bands  in  the  infrared  studied  as  a 
function  of  the  temperature  allows  to  find  the  critical  points  of  the 
combined  density  of  states.  It  should  be  mentionel  that  this  critical 
points  do  not  have  to  correspond  always  to  the  critical  points  for  the 
density  of  states  for  simple  phonons.  The  condition  grad  w(-ij)  =  0 

can  be  satiefyed  when  grad  w(^)  =  -  grad  with  +  for  subtrac¬ 

tive  and  -  for  additive  processes  which  correspond  to  equal  and  op¬ 
posite  of  sign  slops  of  the  two- phonon  branches.  It  is  also  satisfied 
when  grad  w(^)  =  grad  -  0  which  corresponds  to  critical  points 

of  the  density  of  states  curve  for  simple  phonons. 

It  cannot  be  distinguished  between  this  two  cases  experimentally. 
It  is  then  necessary  to  perform  the  calculation  of  the  dispersion 
curves  on  the  basis  of  a  given  model  and  calculate  afterwards  the 
hystogram  of  the  dbuble-phonon  density  of  states  and  compare  this 
hystogram  to  the  ers>erimental  curve.  In  the  case  of  a  rigid  ion  model 
for  CdS  this  work  has  been  done  recently  by  Nusimovici  and  Birmon  . 

The  qualitative  agreement  between  the  two  curves  is  satisfactory.  This 
analysis  shows  also  that  all  the  early  works  on  multiphonon  spectra 
implies  a  great  part  of  arbitrary  assignments. 


J 
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In  Pig. 17  are  presented  on  the  same  fraquenoy  scale  the  two- 
phonon  density  of  states  and  the  absorption  spectrum  corresponding 
to  double-phonon  transitions  in  CdS. 


Comparison  between  ths  calculated  and  tha  experimental 
double  phonon  densities. 
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5.  Localized  Modes  of  Vibrations. 

5.1.  Localized  Surface  Modes. 

For  a  finite  crystal  the  surface  atoms  vibrates  in  a  manner 

14 

different  from  the  atoms  in  the  volume.  Wallis  has  studied  di¬ 
atomic  cubic  crystals  and  find  that  in  addition  to  the  normal  modes 
of  vibrations  it  appears  localized  surface  modes  whose  frequencies 
are  situated  in  the  forbidden  frequency  gap  between  the  acoustical 
and  optical  bands.  The  surface  modes  have  non-vanishing  electric  di¬ 
pole  moments  and  may  give  rise  to  infrared  absorption.  Such  absorp¬ 
tion  has  not  yet  bean  observed  experimentally. 

In  order  to  be  able  to  observe  localized  surface  modes  it  is 

necessary  that  the  ratio  surface/volume  is  as  large  as  possible.  If 

one  considers  thin  film  samples  the  surface  atoms  can  be  seen  as  in- 

14 

dependent  impurities  atoms.  There  are  approximatively  10  atoms  per 

p 

cm  at  the  surface  and  for  a  layer  of  thickness  d  this  would  ccrre- 
14  /  5 

spand  to  10  /d  per  cm  impurity  atoms.  Localized  modes  for  impur- 

18  -3 

ities  become  observable  for  concentration  of  the  order  of  10  cm 
and  this  means  that  surface  modes  can  be  observed  in  thin  films  of 
thickness  1  |i  . 

5.2.  Localized  modes  of  vibrations  due  to  impurity  atoms  in  II- VI 
compounds. 

Localized  modes  have  been  observed  by  means  of  infrared  spectro- 

15 

scopy  in  homopolar  semiconductors  like  Si  containing  light  impurities  , 

1 6  17 

in  intermetallic  semiconductors  like  GaAs  with  Li  and  P  *  and  inSb 

with  Al  as  well  as  in  ionic  crystals  The  frequencies  of  the 

localized  modes  due  to  substitutional  atoms  and  their  coupling  with 

20 

the  radiation  field  can  be  deduced  on  the  basis  of  simple  theory 
when  no  changes  in  force  constant  are  involved. 

Some  semiconductor  crystals  of  the  II-VI  group  such  as  ZnSe  and 

21 

CdS  doped  with  Li  and  Al  have  been  investigated  reoently.  Infra¬ 
red  absorption  due  to  localized  vibrations  in  CdSe  containing  S  and 

22 

CdTe  containing  Se  and  Li  has  also  bsen  investigated  .  Here  the 
frequencies  of  the  absorption  peaks  attributed  to  localized  modes 

23 

are  compared  with  that  calculated  for  a  mass  defect  in  wurtzite  CdS 
subject  to  appropriate  normal  modes  frequency  scaling. 


*****  .  WMtaMMMfcjNi  ^ 
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5.2.1.  Theoretical  considerations. 

The  perfect  lattice  vibrational  frequencies  o>  are  roots  of  the 
secular  equation 

Det  |Muj2-«|  =0  (5.1) 

where  $  is  the  potential  energy  matrix  for  the  perfect  crystal  formed 
by  atoms  with  mass  M.  The  introduction  of  an  imperfection  into  the 
crystal  brings  a  perturbation  c  to  the  potential  energy.  The  matrix 
representing  the  potential  energy  of  vibration  of  the  crystal  with 
impurities  is  *+c.  Tbe  secular  equation  for  the  perturbed  crystal  is 


Det  I  Mu -«-c I  =  0 


(5.2) 


It  is  possible  to  define  a  matrix  g  whose  coefficients  are  the 
Green '8  functions  of  the  crystal  : 


g(Mtd -<&-c)  =  1  -  gc  . 


(5.3) 


The  localized  modes,  due  to  presence  of  impurities,  are  then  solu¬ 
tions  of  the  secular  equation 


Det  1 1  -  g(u>)  c(w)|  =  0 


(5.4) 


In  the  wurtzite  structure  CdS  the  unit  cell  1  contains  four 
atoms  k  =  1,  2,  3,  4  which  are  respectively  Cd^,  Cd^,  Sj,  Sjj  . 

We  just  consider  a  single  mass  defect  (the  force  constants  are  left 
unchanged).  The  defect  subspace  is  three  dimensional  and  localized 
mode  frequency  is  a  solution  of  the  secular  equation  (5.4). 

The  matrix  c  is  a  unit  three  dimensional  matrix  multiplied  by 
(M^-M' )ur  where  M'  is  the  impurity  mass  and  the  mass  of  the  re¬ 

placed  host  ion. 


The  matrix  g  Green's  functions  elements  are  deduced  from  the 


phonon  eigenvectors  j(ja)  and  eigenfrequencies  ur(ja)  s 

%nV  a i V i ( ^  *  ^  )(u)-2_  4 /?  r  2  2/  \i 

*  ig  N  “k*  Cto2-aj2(3a)] 


(5.5) 
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Prom  the  point  group  symmetry  of  the  crystal  around  the  defect  the 
Green's  functions  satisfy  the  following  relations  : 


gak,a'k'(°’0)  “  6oa'  gak,ak^0’0^ 


(5.6) 


gxktXk(°,0)( «2)  =  gyk,yk(°.0)  (a»2) 


The  impurity  may  replace  either  Cd  or  S  and  four  different  Green's 
functions  are  used  in  this  problem  : 


«.  -  “l  621,2,(o.=>(“2>  . 


(5.7) 


*3  '  "3  ®Z3,Z3(0,0)(“  >  > 


*♦  ‘  “3  **3,r3(0-0)("  > 


where  M1  and  are  respectively  the  mass  of  Cd  and  the  mass  of  S. 


For  many  defect  vibration  problems  and  especially  for  band 
resonance  calculations  it  is  convenient  to  know 


lim  <p  -  0  =  g'(w)  +  2S 


(5.8) 


The  values  of 


vak,a'k'*1,1,^w*  “  n  H  aak^a^  Vk'^a)  exp  ia(rl-rl')  &[<*>(&)-<*>] 

da 

(5.9) 

are  obtained  by  using  the  computed  CdS  spectrum  *  for  125  values  of 
a  in  the  reduced  Brillouin  zone.  The  frequency  band  (w  \  u^j)  is  di¬ 
vided  into  Z  =  61  equally  spaced  regions  and  the  sum  (5.9)  carried 
over  m(ui^/Z)  K.  u>  (  (m+l)(u^Z)  to  from  a  hystogram  for  v.  The  va¬ 
lues  of  g'(w)  are  then  obtained  at  the  mid  points  of  the  hystogram 
sections  by  using  the  Kramer s-Krbnig  relation 

rt‘l®  Vj  (u' )  dw' 

«*<«>-]  —z — 1  • 

1  {  ur  -  a'* 


as 


Flgura  20.  Rati  part  of  tha  Graen’a  functions. 

2  2 
w  g'<(u)  correaponda  to  curve  1,  u  g'  («)  to  curve  2, 
2  A  2  2 

*>  g'^Cw)  to  curva  3,  and  u  g’^  (u)  to  curva  4. 


The  localized  mode  frequencies  are  solutions  of  the  following 
equations  : 


e1  a)  g1  =  i  for  the  mode 

el  gg  =  1  for  the  mode 

where  the  impurity  replaces  Cd  and 
2 

ui  =  1  for  the  mode 
2 

Ej  u  g^  =  1  for  the  r5  mode 
where  the  impurity  replaces  S 

e1  =  (M1-M')/M1  and  =  (M3-M')/M3 


(5.10) 


(5.1l) 


--*•  *  •  - j_i v/  L*m,e**\  • 


■j-^awi.r,ni 
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The  imaginary  part  of  the  Green's  functions  v.M  may  be  com¬ 
pared  with  the  one-phonon  density  of  states  in  CdS  recently  calculated 
and  shown  in  Fig.21.  The  functions  v^u)  and  v2(u)  related  to  the 
heavy  Cd  ion  displacement  vectors  are  enhanced  in  the  low  frequency 
acoustic  region  and  the  functions  v^u)  and  v4(«)  related  to  the  light 

S  ion  displacement  vectors  are  enhanced  in  the  high  frequency  optic 
region. 


L*g.ure  ,?1»  One  phonon  density  of  states  for  CdS 
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a)  Localized  Si  ode  a. 

The  localized  node  frequencies  obtained  from  equations  ( 5  - 10) 
and  (5.1l)  are  pi  at  ted.  in.  Pig. 22  versus  e. 


Figure  22.  Localized  mode  frequencies  for  a  mass  defect  — -  -  1 

is  plotted  versus  for  curves  1  and  2  where  tfie 
impurity  replaces  Cd,  and  versus  for  curves  3  and  4 
where  the  impurity  replaces  S.  Curves  1  and  3  correspond 
to  rx  modes  and  curves  2  and  4  to  modes. 

If  we  suppose  <1^  =  306  cm  1  for  CdS  a  localized  mode  appears  only  if 
^  O.36  for  an  impurity  replacing  S  and  e1  0.58  far  an  impurity 
replacing  Cd.  If  e  is  large  enough,  two  modes  appear  out  of  the  band; 
a  singlet  and  a  doublet  with  a  higher  frequency  than  the  singlet. 
The  splitting  increases  when  the  impurity  becomes  lighter.  Up  to  now 

localized  modes  due  to  a  light  impurity  mass  defect  have  only  been 

c.  n  pa 

observed  for  Liw  and  Li  in  CdS  .  We  shal  discuss  such  modes  observ¬ 
ed  in  other  II-VI  semiconductors.  Predictions  for  different  impurities 
in  CdS' are  summarized  in  Table'VI. 

Table  VI 

Localized  mode  frequencies  for  substitutional 
impurities  in  CdS 


Impurity 

Replaced  Ion 

M1 

£ 

o>ri  cm-1 

wIV  cm" 

5 

Be 

Cd 

9 

0.92 

540 

580 

Mg 

Cd 

24 

0.786 

365 

380 

Ca 

Cd 

40 

0.683 

320 

330 

0 

S 

16 

0.5 

324 

339 

Li 

Cd 

7 

0.94 

615 

665 

i 


For  light  impurities  (e  -  l)  all  the  localized  mode  frequencies 
tend  to  infinity.  This  can  he  shown  from  the  limit  form  of  the  Green's 
functions  vhen  oi  is  much  bigger  than  (11(33.)  : 

2 

,  ..  o  _ _  I  3si)  I  . 

gok,ok( )  ~  II  5  —“2  * 

ja  \(i) 

but 

0ok^a^  apk'^a^  =  6«fl  » 


aP  kk' 


then  when  u>  -»  00 


and  the  solution  for  e  is  e  =  1  . 


b)  Band  Modes. 

When  a  heavy  mass  impurity  is  introduced  a  band  resonance  may 
appear  at  oi  \  The  response  of  the  lattice  is  related  to  the  imag¬ 
inary  part  of  the  imperfect  crystal  Green's  functions 

[Kv,  ((i))]/2w 

Im  g,  (di)  =  - - - 5 - - -  . 

1  [l  -  em2g£(a>)]2  +  [TIeuvi(ai)]2/2 

So  a  resonance  is  well  defined  at  the  <i>  which  is  a  solution  of 

1-eoi  g'(oi)  =  0  only  if  v,  (oi)  is  small.  For  3e  in  CdS,  e  =  -1.5  and 

1  1  - 1 
a  resonance  is  expected  at  210  cm 


The  model  on  which  this  calculation  is  based  being  simple,  one 
is  justified  to  extend  the  results  to  other  II-VI  compounds  when  no 
considerable  changes  in  farce  constants  are  expected  when  going  from 
one  compound  to  another. 


In  the  cubic  crystals  of  this  group,  like  CdTe  for  ex&mple,  the 
site  symmetry  is  T^  and  therefore  one  triply  degenerate  mode  is  to 
be  expected. 


then  introduced  in  CdS  single  crystals  the  Se  with  mass  M  =  79 
replaces  the  S  atom  with  mass  m^  =  32,  Cd  having  a  mass  of  m1  =  112.4, 
we  have  the  relation  ^  H  K  .  The  mass  parameter  in  this  case  is 
e  *  (mrM)/m1  =  -1  .47  .  The  infrared  reflectivity  spectrum  (Pig. 23)  of 
this  system  shows  a  resonance  at  ISO  cm-1. 


Figure  23. 


Infrared  reflectivity  epectrua  of  CdS  with  1.5  X  Se 
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This  frequency  being  out  of  the  rcstraMen  region  the  corresponding 
mode  could  be  considered  as  a  gap  mode.  The  density  of  mode  at  this 
frequency  (see  Pig. 21 )  is  very  low,  practically  zero,  one  can  there¬ 
fore  consider  that  this  mode  is  situated  in  a  forbidden  frequency 
region  ,  the  frequency  gap  ,  and  is  therefore  called  a  gap  mode. 

With  increasing  concentration  of  Se  the  gap  mode  frequency,  <n  =  180 


-1 


g 


cm  ,  splits  off  into  two  modes  having  respectively  the  characters 
of  a  transverse  optical  and  a  longitudinal  optical  modes. The  frequen¬ 
cies  of  these  two  modes  fall  apart  from  each  other  linearly  with  the 


concentration  and  reach  the  limit  of  and  for  CdSe  respectively 
=  170  cm-1  and  0Lq  =  210  cm-1 


TO 


when  Se  approaches  50^.  The  varia¬ 
tions  of  and  (jTq  deduced  from  the  position  of  each  concentration 
of  crossing  of  the  zero  line  for  the  real  part  of  the  dielectric  con¬ 
stant  in  changing  from  negative  to  positive  values  and  from  the  posi¬ 
tion  of  the  peak  of  the  imaginary  part  of  the  dielectric  constant  are 
shown  in  Pig. 24. 


Figure  24.  Frequency  variation  of  TO  and  L0  in  the  CdS-Se  alloy 

with  increasing  concentration  of  Se . 

0  A 

In  earlier  investigations  ^  concerning  the  reflectivity 
spectra  of  CdS-Se  alloys  analogous  results  have  been  obtained. 


39 


5. 2. 2. 2.  S  in  CdSe. 

For  concentration  higher  than  50$  of  Se  in  CAS  the  situation 
will  be  inverted,  the  inpurity  will  be  now  S  in  CdSe.  S  is  a  light 
impurity  M  =  32  in  CdSe  substituting  Se  with  m^  =  79 >  the  defect  mass 
parameter  e  here  is  E  =  (79-32)/79  =  0.595. 

In  the  wurtzite  structure  CdS  have  a  unit  ceil  1  containing 
four  atoms  k  »  1,  2,  3»  4,  which  are  respectively  Cd^,  Cd^,  Sj  and 
Sjj.  The  point  defect  introduced  in  the  perfect  crystal  substitutes 
one  of  these  atoms.  In  a  first  approximation  only  the  mass  defect  is 
considered  assuming  that  substitution  of  an  atom  for  another  does  not 
introduce  changes  in  the  force  constants  of  the  perfect  crystal.  The 
localized  mode  frequency  is  a  solution  of  the  secular  equation  (5.4). 

To  solve  this  equation  it  is  necessary  to  take  only  a  3*3  submatrix 
G(w)  corresponding  to  the  defect  subspace.  Substitution  of  the  Cd 
and  S  atoms  are  studied  with  the  mass  coefficients  respectively 

e1  =  AM/M1  ,  Cj  =  AM/Mj 

In  the  wurtzite  crystal  the  point  defect  has  no  more  the  site  sym¬ 
metry  Tp  but  Cgv  and  the  triply  degenerate  localized  mode  is  now 
split  into  two  modes  1’^  and  simply  and  doubly  degenerate.  The  local¬ 
ized  mode  frequencies  are  solutions  of  the  following  equations  when 

the  impurity  replaces  Cd  atoms  : 

2 

.  e^u>  g^  =  1  for  the  mode  doubly  degenerate  corresponding  to 
vibrations  of  the  impurity  atoms  perpendicular  to  the  c  axis  of  the 
crystal  ; 

2 

.  e^w  g1  =  1  for  the  T1  mode  non-degenerate  corresponding  to 
vibrations  of  the  impurity  atoms  parallel  to  the  c  axis  of  the  crystal. 

When  the  impurity  replaces  the  S  atoms  the  localized  mode  fre¬ 
quencies  are  solutions  of  the  following  equations  : 

2 

EjUrg^  =1  for  the  mode  , 

2 

e^o)  gj  =  1  for  the  mode. 

The  position  of  the  localized  modes  is  obtained  directly  accord- 

p  p 

ing  to  these  equations  from  the  plot  of  oj  g(io  )  as  a  function  of  the 

p  p  | 

frequency  in  taking  w  g(w  )  =  -  . 
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If  we  admit  that  the  density  of  states  does  not  considerably 
changes  from  CdS  to  CdSe  we  may  then  use  the  calculated  curve  for  CdS 
given  in  Pig. 22,  provide  we  make  the  appropriate  frequency  scaling 
in  substituting  for  the  value  of  for  CdSe  which  is  u^q  =  210 


cm 


-1 


to  deduce 


“L 


The  values  obtained  in  this  way  are  co,  ( r  )  =  244 
—  1  —  1  ■*-  * 
cm  and  i^(r^)  =  256  cm  .  The  infrared  reflectivity  spectrum  for 

CdSe  containing  S  is  given  in  Fig. 25. 


FlRure  25.  Infrared  reflectivity  of  CdSe  with  0.5  XS» 

The  experimental  values  deduced  from  Fig. 25  are  respectively 

u)(r,)=  266.5  cm-1  and  o>( r_ )  =  269  cm-1. 

1  5 

The  difference  between  the  calculated  and  observed  values  could 
be  attributed  to  a  decrease  of  the  first  neighbours  force  constants 
when  S  is  substituted  for  Se. 

The  line  width  is  of  the  order  of  4  cm-1  and  the  relative  line 
width  hu/u  =  0.015  which  is  of  the  order  of  magnitude  for  relative 
broadening  when  Au>  describes  the  broadening  of  an  optical  spectral 
line  due  to  enharmonic  interaction  involving  the  vibrational  energy  of 
the  local  mode.  A  temperature  line  shift  is  observed  of  the  order  of 
3  cm-1  between  25 °K  and  300°K. 

On  the  basis  of  symmetry  consideration  it  is  easy  to  deduce  the 
selection  rules  showing  that  from  the  split  of  and  r modes  only 
one  will  be  optically  active  for  polarized  light.  The  r  will  be  active 
for  light  with  E//C  and  will  be  active  for  eJ_C.  It  is  surprising  to 
find  out  that  experimentally  the  two  picks  seem  to  be  optically  active. 


1 
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When  the  concentration  of  S  in  CdSe  increases  the  two  localized 
modes  shifts  in  frequency  until  reaching  the  limit  frequencies  corre¬ 
sponding  to  the  normal  modes  in  CdS  respectively  uj^q  =  240  cm  1  and 
~  300  cm  1 . 

6. 2. 2. 3.  light  impurities  in  CdTe. 

CdTe  has  the  cubic  zinc-blende  (sphalerite)  structure  band  on 
2 

the  cubic  space  group  Tt.  The  dispersion  curves  calculated  for  a 
hypothetical  cubic  CdS  show  only  slight  difference  from  that  of  wurtz- 
ite.  If  it  is  acceptable  to  use  the  same  frequency  distribution  as  in 
CdS  it  will  be  possible  with  an  appropriate  frequency  scaling  for  LO 
to  deduce  a  value  for  for  each  impurity.  The  frequency  of  the 
localized  mode  obtained  in  this  way  may  represent  a  good  order  of 
magnitude  to  compare  with  the  exprimental  results. 

a)  Se  in  CdTe. 

Se  M  -  79  is  a  light  impurity  in  CdTe  which  substitutes  to  the 

Te  atom  whose  mass  is  m  =  127.6.  The  mass  parameter  in  this  caso  is 

^  -  1 
E  =  0.38.  The  longitudinal  optical  mode  frequency  is  =  16?  cm  . 

Prom  the  relations  given  in  Pig. 22  it  is  easy  to  deduce  the  value  of 

the  localized  mode  frequency  which  is  (^  =  167  cm-1. 

The  infrared  reflectivity  spectrum  for  CdTe  containing  Se  is 
given  in  Pig. 26. 


Figure  26 


Infrared  reflectivity  spectrum  of  CdTe  with  a,  1  I  3« 


-#• 
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The  Kramers-Kronig  analysis  of  the  optical  constants  is  repre¬ 
sented  graphical!,'  in  Fig. 27. 


From  the  peaks  of  the  imaginary  parts  of  the  dielectric  con¬ 
stant  we  deduce  the  values  of  u_n  =  139  cm”1  and  the  localized  mode 
frequency  =  170  cm  .  The  calculated  Im  -r  shows  two  peaks  at  fre¬ 
quencies  162  and  180  cm-1.  For  pure  CdTe  the  longitudinal  optical  mode 
frequency  is  ui^q  =  167  cm_1;in  doped  crystals  we  have  =  162  cm-1 
and  =  188  cm_  1 »  “lq  woult*  then  3Uf f er  a  slight  shift  of  5  cm” 1 

when  a  light  impurity  is  introduced. 


The  line  width  for  is  of  the  order  of  10  cm-1  and  the  re¬ 
lative  line  width  du/w  -  0.059  which  can  still  be  considered  as 
relative  enharmonic  broadening. 


b)  Li  in  CdTe. 


The  introduction  of  Li  in  CdTe  has  been  done  by  the  method  of 
progressive  crystallization.  The  crystals  obtained  by  this  method 
when  at  low  concentration  of  Li  approximately  0.1$  are  n-type  and 
with  low  free  carrier  concentration  of  the  order  of  SxiO1^  cm"'*.  At 
higher  Li  concentration  approximately  1$  the  samples  are  of  p-typef 
the  free  carrier  concentration  is  very  low  Qf  the  order  of  5x  1 0 1  ^cm" •* . 
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The  samples  with  low  Li  concentration  are  opaque  in  the  infra¬ 
red.  At  higher  Li  concentration  the  samples  become  transparent  proba¬ 
bly  because  of  a  partial  compensation. 

At  low  concentration  Li  is  supposed  to  occupy  interstitial 
sites  and  tc  act  as  a  donor  whereas  at  high  concentration  the  occupa¬ 
tion  of  substitutional  sites  where  Li  acts  as  an  acceptor  becomes 
predominant.  In  the  neighbourhood  of  1 $  of  Li  a  partial  compensation 
is  achieved  and  the  samples  are  transparent.  The  absorption  spectra 
of  crystals  with  of  Li  at  300°K  and  100°K  are  given  in  Pig. 28. 


Figure  28.  Absorption  spectra  of  CdTe  with  1  %  Li  at  300*K 

Pure  crystal  - 

1  *  LI  —  .  —  .  — 

and  at  100*K 

Pure  crystal  ——————— 

1  X  LI  —  ..  —  ..  — 

These  spectra  clearly  show  absorption  peaks  due  to  the  presence  of  Li. 
The  largest  peak  is  centered  at  270  cm-1  at  300°K  and  it  shifts  to¬ 
ward  the  high  frequency  for  lower  temperatures  :  the  maximum  of  the 
peak  is  situated  at  275  cm-1  when  taken  at  100°K  and  it  is  at  277  cm"1 
at  25 °K. 

The  width  6f  this  band  of  the  order  of  35  cm”1  is  anormally 
large  so  is  also  the  relative  band  width  Aw/w  =  0.130. 


i 
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5.2. 2. 4.  Impurities  in  ZnSe. 

ZnSe  is  again  a  cubic  zinc-blende  crystal  for  which  the  analysis 
of  the  infrared  spectra  concerning  the  phonon  interaction  with  the 
radiation  field  can  be  treated  on  the  basis  of  the  equivalent  disper¬ 
sion  relation  for  the  sphalerite  CdS.  The  tacite  assumption  in  calcula¬ 
tion  the  local  mode  frequencies  is  that  the  frequency  distribution 
i3  analogous  to  that  of  CdS.  The  experimental  results  are  very  close 
to  that  calculated  in  this  hypothesis  which  can  be  therefore  consider¬ 
ed  as  a  satisfactory  first  order  approximation. 

a)  Mn  in  ZnSe. 

When  Mn,  M  =  55 t  is  introduced  in  ZnSe  it  substitutes  the  Zn 
atom  of  mass  =  65.4.  The  mass  parameter  is  then  e  =  0.'6.  The  mass 

difference  being  extremely  small  e  K.  0.30  for  a  =  25?  cm”1  the  mode 
due  to  the  presence  of  this  in  purity  is  a  band  mode  and  therefore  not 
directly  observable  in  the  frequency  range  out  of  the  restrahlen. 

The  infrared  reflectivity  spectrum  of  this  system  is  given  in 
Fig.  25. 


Figure  29. 


Infrared  reflectivity  of  ZnSe  with  1.5  X  Mn. 
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A  close  examen  of  the  optical  constants  displayed  in  Pig. 30 

shows  a  structure  on  the  imaginary  part  of  the  dielectric  constant 

TO 


which  gives  two  frequencies  for  the  transverse  optical  mode 

19^  cm-1  and  =  202  cm-1. 

TO 


:  ^ 


Figure  30.  Dielectric  constants  of  ZnSe  with  1.5  X  Mn. 

For  the  pure  ZnSe  we  have  uT0  -  202  cm"1.  It  seems  that  the  presence 
cc  an  impurity  shifts  the  u^q  from  ■  252  cm  1  to  an  higher  energy 

“lo  “  254  an_1 

b)  S  in  ZnSe. 

S  is  a  light  impurity  in  ZnSe  substituting  Se.  The  mass  factor 
is  e  -  0.595.  The  frequency  scaling  done  with  oj^q  =  252  cm”1  leads  to 
a  localized  mode  frequency  =  297  cm  1 . 


The  reflectivity  spectrum  for  3  doped  ZnSe  is  represented  in 


Fig. 31 • 


Figure  31.  Infrared  reflectivity  of  ZnSe  with  »v/ 1  Z  S. 

The  real  and  imaginary  parts  of  the  dielectric  constant  deduced 
hy  the  Kramers-Kronig  analysis  are  desployed  in  Fig. 32. 


Figure  32.  Dielectric  constants  of  ZnSe  with  /v  1  Z  S. 
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normal  mode  frequencies  for  fhe  undoped  crystal  are  =  252  cm 


The 


u>T-  =  202  cm-1 
=  297  cm'1 


,  and  for  the  S  doped  the  local  mode  frequencies  are 

and  idj10**- 305.5  cm"1- 


The  line  width  which  is  function  of  the  concentration  is  in  this 

-1 


case  (~i^  S)  Aw  =  22  cm 
0.074. 


and  the  relative  line  broadening  is  Aw/o>  = 


c)  Te  in  ZnSe . 

Te  with  M  =  127.6  is  a  heavy  impurity  in  ZnSe.  The  mass  para¬ 
meter  is  e  =  -0.615.  A  heavy  defect  in  the  lattice  should  lead  to  a 
resonance  in  the  frequency  allowed  region  or  to  a  gap  noae  if  the 
density  of  states  at  the  given  frequency  is  very  low  or  pratically 
zero. 

The  infrared  reflectivity  spectra  for  pure  and  Te  doped  ZnSe 
are  given  in  Fig. 33. 


Figure  33. 


Infrared  reflectivity  of  ZnSe  with  4  X  Te , 


i 
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The  real  and  imaginary  parts  of  the  dielectric  constant  are 
deaployed  in  Fig. 34. 


Figure  34.  Dielectric  constants  of  ZnSe  with  4  *  Te . 

From  the  difference  between  the  imaginary  part  of  the  dielectric 

constant  for  pure  ZnSe  and  Te  doped  ZnSs  one  can  deduce  the  frequency 

of  the  gap  mode  u>  =  182  cm”1.  The  width  of  this  band  is  Au  =  40  cm”1 
€ 

which  combines  the  impurity  line  width  and  the  band  broadening  due  to 
the  concentration.  The  relative  band  broadening  is  Aw/w  -  0.22  . 
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7.  Discussion. 


In  spite  of  the  considerable  amount  of  work  presently  under  way 
in  the  lattice  dynamics  of  the  II-VI  compounds  and  the  considerable 
progress  it  still  remains  much  to  be  done. 

The  dispersion  relations  for  OdS  have  teen  calculated  in  terms 
cf  an  extremely  simplified  model  whose  physical  significance  remains 
to  be  explicited.  Although  the  first  results  are  satisfactory  and  the 
eigenvalues  and  eigenfunctions  have  served  as  basis  for  the  calculation 
of  the  local  mode  frequencies  which  are  in  good  agreement  with  the 
experiment,  it  is  certainly  suitable  t'  improve  this  model  including 
physical  parameter  such  a3  the  electron  polarizability  of  the  ions. 

The  calculated  double-phonon  density  of  states  do  not  allow  a  direct 
comaprison  with  the  two-phonon  optical  spectra  in  the  whole 
spectral  range.  Such  a  descrepenoy  should  be  understood  eather  in 
terms  on  unsatisfactory  physical  implications  in  the  calculation  of 
the  frequency  distribution  curve  or  in  terms  of  second  order  matrix 
elements  for  optical  transitions  employing  more  than  two  particles. 

The  qualitative  features  of  the  normal  mode  spectra  are  well 
understood  and  the  spectral  measurements  give  satisfactory  quantita¬ 
tive  data  for  those  modes  which  are  coupled  to  the  radiation  field 
at  the  appropriate  points  of  the  Briilouin  zone  which  are  eather  the 
center  or  high  symmetry  edges. 

Studies  on  thin  films  have  proved  to  be  extremely  useful  for 
the  experimental  determination  of  the  normal  mode  frequencies. 

Surface  modes  are  still  to  be  seen  experimentally. 

Elementary  calculations  using  Green  function  technics  and 
assuming  isotopic  substitution  without  force  constant  changes  have 
proved  to  be  an  extremely  useful  guide  in  the  discussion  of  the 
experimentally  determined  localized  modes  due  to  impurities.  The 
orders  of  magnitude  calculated  even  in  drastic  assumptions  coincide 
very  closely  to  with  the  experimentally  measured  values.  The  reason 
is  probably  the  fact  that  the  calculated  density  of  states  obtained 
in  a  very  rough  model  which  is  used  for  the  calculation  of  the  local¬ 
ized  modes  frequencies  is  insensitive  to  3light  changes  in  force  con- 
tants  when  an  ion  is  replaced  with  a  similar  ion  from  the  same  groups. 
Only  an  improved  initial  model  may  lead  to  results  more  sensitive  on 
subtile  force  constant  changes  produced  by  the  substitution  of  analog¬ 
ous  ions. 
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The  first  measured  localized  mode  frequencies  in  II-VI  compounds 

pi 

are  those  concerning  the  Li  doped  ZnSe  and  CdS  .  In  the  case  of  CdS 

doped  with  Li  anormally  large  bands  have  been  observed  and  attributed 

7  _ «  6  _  i 

to  the  presence  of  Li  at  457  cm  and  Li  at  474  cm 

If  Li  wa3  substituting  Cd  atom  the  mass  difference  parameter 
7  6 

would  be  e  =  0.938  for  Li'  and  t  =  0.946  for  Li  ,  this  would  give  in 
tha  simplest  approximation  of  iootcpic  mass  substitution  for  the 
localized  mode  frequency  of  the  impurity  atom  vibration  along  the  c 
axis  the  following  values  : 

7  - 1 

for  Li  ,  u>p1  =  492  cm  , 

f  or  Li^  ,  wp 1  =  503  cm- 1  . 

The  calculated  values  are  larger  than  that  observed  experiment¬ 
ally  there  would  indicate  that  Li  interacts  stronger  with  the  host 
lattice  than  a  simple  mass  defect.  A  change  of  force  constant  would 
explain  this  difference.  Increase  in  the  deviation  of  force  constant 
from  the  host  lattice  force  constant  acts  in  the  same  direction  as  a 
reduction  of  mass  i.e.  increase  of  the  localized  mode  frequency. 

The  presence  of  3  in  CdSe  crystals  introduces  a  mass  defect 
parameter  e  =  0.595  and  hence 

Wp1  =  244  cm-1  ,  =  256  cm-1 

Experimentally  two  peaks  are  observed  at  266.5cm-1  and  269  cm  1  at 
room  temperature.  The  calculated  values  now  are  smaller  than  the  ob¬ 
served  ones.  This  discrepency  could  be  also  attributed  to  a  change 
in  force  constant  where  Se  is  replaced  by  S,in  this  case  the  substi¬ 
tution  would  means  a  decrease  of  force  constant. 

The  uine  width  is  of  the  order  of  4  cm-1  and  the  relative  lins 
width  Aio/o)  is  O.O15  which  is  of  the  order  of  magnitude  for  relativs 
broadening  of  the  defect  mode  frequency  when  Au>  describes  the 
broadening  due  to  anharmonic  interaction  of  an  optical  spectral  line 
involving  the  vibrational  energy  of  the  local  mode. 

A  temperature  line  shift  is  observed  of  the  order  of  3  cm-1 
between  25 °K  and  309°K. 

In  the  cass  of  Li  doped  CdTe  an  absorption  band  is  observed 
centered  at  270  cm-1  which  could  be  attributed  to  localized  modes  due 
to  ths  presence  of  Li  but  whose  interpretation  is  difficult  for  many 
reasons  a  i  particularly  because  of  the  band  width  which  is  of  the 
order  of  35  cm-1  and  does  r.ot  seem  to  vary  considerably  with  temper¬ 
ature.  The  relative  band  width  Aw/w  here  is  0.1 30  what  is  by  a  factor 
of  10  larger  than  what  would  be  expected  in  the  case  of  single  subs¬ 
titutional  impurity  vibrations. 
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The  essential  contribution  to  the  residual  line  width  when 
T  **  0  comes  from  the  contribution  of  the  spontaneous  decay  of  the 
excited  localized  vibrations  into  band  modes.  In  this  process  the 
energy  conservation  principal  XI  ^  =  u  should  he  respected.  In  the 
case  of  Li  in  CdTe,  as  well  as  in  the  other  cases  we  have  examined 
here,  u  K.  2  ,  consequently  the  decay  of  a  localized  mode  into  two- 

band  modes,  is  possible  through  the  third  order  anharmonic  potential 
and  this  decay  could  be  the  dominant  process  at  room  temperature2^’2^. 

In  order  to  understand  the  possibility  to  obtain  optically 
transparent  CdTe  with  relatively  high  concentration  of  Li  we  had  to 
suppose  that  at  low  concentration  Li  is  donor. and  at  high  concentra¬ 
tion  gives  rise  to  acceptor  levels.  This  leads  to  the  possibility  of 
an  almost  complete  compensation  for  a  given  concentration  of  dopant. 

The  interaction  of  an  interstitial  atom  with  the  lattice  is  not  the 
same  as  that  on  a  substitutional  impurity.  It  is  therefore  clear 
that  although  the  mass  is  the  same  due  to  changes  inthe  force  con¬ 
stant  the  frequencies  of  the  localized  modes  should  be  different  for 

the  two  situations.  One  should  in  such  a  case  observe  two  absorption 

6  7 

peaks.  In  addition  Li  has  two  isotopes  Li  and  Lx.'  constituting  each 
one  a  mass  defect  and  therefore  giving  rise  to  an  absorption  band  each. 
All  these  are  reasons  to  believe  that  the  large  band  at  270  cm-1 
contains  a  certain  number  of  peaks  remaining  to  be  resolved  in  the 
simplest  case. 

Experiments  are  now  undertaken  at  different  concentrations  of 
Li  in  order  to  further  investigate  this  situation. 

The  presence  of  Se  in  CdTe  leads  to  a  typical  localized  mode 
whose  frequency  deduced  from  the  Kramers-Kr'dnig  analysis  of  the  re¬ 
flectivity  spectrum  in  ux.  =  170  cm-1  whereas  the  calculated  value 

i 

with  the  appropriate  normal  mode  frequency  scaling  is  =  i 67  cm  . 
The  presence  of  impurities  at  relatively  high  concentration  seems  to 
have  an  effect  on  the  longitudinal  mode  frequency  which  is  expressed 
by  a  slight  shift  of  the  normal  mode  freqxiency  is  now  162  crn” 
and  a  longitudinal  part  to  the  local  mode  frequency  appears  at  180 
cm-1.  At  high  impurity  concentration  one  can  therefore  define  a  trans- 
verse  and  a  longitudinal  part  for  the  local  mode  ux^  =  l7Q  cm  and 
ox.  x,°  =  180  cm-1  che  calculated  value  for  the  local  mode  frequency  is 
within  these  limits  :  ux^  =  167  cm 

Analogous  situation  has  been  formed  in  ZnSe  doped  with  S  where 
the  calculated  value  ox^  =  297  cm-1  is  within  the  limits  of  that  ob¬ 
tained  experimentally  :  =  297  cm  1  and  ux^^  *  3Gv5om  1 . 
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Gap  modes  or  band  resonances  have  also  been  observed  in  the 

case  of  heavy  impurities  as  for  example  ?e  in  ZnSe  where  u  =  182  cm 

£ 

has  been  deduced  from  the  analysis  of  the  reflectivity  spectra. 


A  quite  complete  investigation  on  the  effects  of  heavy  doping 

has  been  carried  in  the  case  of  the  system  CdSe-S.  For  Cd3  containing 

low  concentration  of  Se  the  dopant  acts  as  a  heavy  impurity  and  gives 

a  gap  mode  at  <d  =  ISO  cm-1.  The  increase  of  the  3e  concentration 
& 

brings  up  two  constituants  of  the  o>  band  having  respectively  the 

LO  ®  TO 

character  of  a  longitudinal  w  and  a  transverse  mode  u>  .  The  fre- 
LO  TO  S  6 

quencies  w  and  w  are  linearly  dependent  on  the  impurity  concen- 

tration  and  tends  toward  the  normal  modes  frequencies  of  CdSe  at  tne 

high  concentration  limit. 


S  in  CdSe  gives  two  localized  modes  <m.(r  )«  266.5  cm-1  and 
u^dV)  =  269  cm  whose  frequencies  are  also  concentrations  dependent 
and  reach  at  high  S  concentration  the  normal  modes  frequencies  for  CdS. 


It  seems  experimentally  that  we  have  a  continuous  relation 
between  localized  and  normal  modes  frequencies  depending  on  concentra¬ 
tion.  It  will  be  highly  interesting  to  undersee  theore'-*  oally  such  a 
continuous  relation  when  we  pass  from  a  perfect  lattice  to  a  desorder- 
ed  system. 


8.  Conclusion. 

The  discussion  of  lattice  dynamics  of  II-VI  compounds  is  quite 
general  and  interesting  in  one  of  the  fact  that  the  bonding  spanes 
the  whole  range  from  covalent  to  icnic.  A  complete  calculation  of  the 
dispersion  relations  is  now  available  based  on  a  simple  model  in  which 
short  range  forces  simulate  covalent  band  contribution  and  Coulomb 
forces  are  functions  whose  origin  is  the  ionic  nature  of  the  bounding. 
The  parameters  used  in  the  calculation  are  plotted  with  the  optically 
measured  normal  modes  frequencies  at  the  center  of  the  Brillouin  zone. 
The  calculated  phonon  frequencies  of  the  edge  of  the  Brillouin  zone 
in  the  principal  symmetry  directions  compares  well  with  the  experiment¬ 
al  values. 


.  ■  item* 
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Analysis  of  the  coupling  of  the  radiation  field  with  the  phonon 
field  is  developed  in  the  frame  of  a  macroscopic  theory  which  account 
well  for  all  the  experimental  situations.  Measurements  in  blende  and 
wurtzite  structures  are  consistently  analysed. 

It  is  shown  that  numerous  arbitrary  assignments  in  the  multi¬ 
phonon  optical  spectra  can  be  easily  erroneous  because  while  it  is 
clear  that  critical  points  in  the  combined  density  of  states  give  rise 
to  absorption  band  it  is  not  necessary  that  this  critictJ.  point  cor¬ 
responds  to  critical  point  in  the  single  phonon  density  of  states.  In 
order  to  have  a  critical  point  in  the  combined  two-phonon  density  of 
stetes  the  condition  grad  w(q)  =  0  can  be  satisfied  also  when 
grad  w(^)  =  ±  grad  w(^,)  which  corresponds  to  equal  and  opposite  of 
sign  slops  in  the  two-phonon  branches.  It  is  then  cl' ar  that  the  only 
correct  assignments  will  come  from  the  comparison  of  the  optical  ab¬ 
sorption  spectra  with  the  calculated  density  of  states  when  the  dis¬ 
persion  relations  are  known.  An  example  concerning  CdS  is  given  here. 

The  introduction  of  imperfections  in  the  crystal  depending  on 
the  impurity  mass  gives  rise  to  localized  modes,  to  resonances  or  to 
gap  modes.  The  examen  of  a  certain  number  of  systems  has  shown  that 
at  high  impurity  concentration  the  localized  as  well  as  the  gap  inodes 
give  rise  to  a  reertraHen  band  in  the  reflectivity  spectra  which  sub¬ 
mitted  to  a  Kramers-Krbnig  shows  contributions  having  transverse  op¬ 
tical  and  longitudinal  optical  characters.  With  the  increase  of  con¬ 
centration  and  forming  the  ternaryalloys  the  frequencies  of  these  two 
parts  depart  from  each  other  and  at  about  50$  in  the  concentration 
reach  the  values  of  the  normal  modes  frequencies  of  the  system.  This 
frequency  shift  seems  to  be  linear  with  concentration  and  offers  an 
interesting  feature  in  the  study  of  the  transition  from  perfect  crystal 
with  low  impurity  ccnoentratiai  to  the  disordered  ternary  system. 

With  the  development  of  the  investigations  of  the  dispersion 
relations  of  perfect  crystals  we  hope  to  see  also  further  studies  con¬ 
cerning  the  disordered  alloys. 
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